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We analyze heavy quark free energies in 2-flavor QCD at finite temperature and the corresponding 
heavy quark potential at zero temperature. Static quark anti-quark sources in color singlet, octet and 
color averaged channels are used to probe thermal modifications of the medium. The temperature 
dependence of the running coupling, aqq{r,T), is analyzed at short and large distances and is 
compared to zero temperature as well as quenched calculations. In parts we also compare our 
results to recent findings in 3-flavor QCD. We find that the characteristic length scale below which 
the running coupling shows almost no temperature dependence is almost twice as large as the 
Debye screening radius. Our analysis supports recent findings which suggest that Xc and tpt are 
suppressed already at the (pseudo-) critical temperature and thus give a probe for quark gluon 
plasma production in heavy ion collision experiments, while J/ip may survive the transition and will 
dissolve at higher temperatures. 
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I. INTRODUCTION 

The study of the fundamental forces between quarks 
and gluons is an essential key to the understanding of 
QCD and the occurrence of different phases which are 
expected to show up when going from low to high tem¬ 
peratures (T) and/or baryon number densities. For in¬ 
stance, at small or vanishing temperatures quarks and 
gluons get confined by the strong force while at high tem¬ 
peratures asymptotic freedom suggests a quite different 
QCD medium consisting of rather weakly coupled quarks 
and gluons, the so-called quark gluon plasma (QGP) 0. 
On quite general grounds it is therefore expected that 
the interactions get modified by temperature. For the 
analysis of these modifications of the strong forces the 
change in free energy due to the presence of a static quark 
anti-quark pair separated by a distance r in a QCD-like 
thermal heat bath has often been used since the early 
work HQ In fact, the static quark anti-quark free en¬ 
ergy which is obtained from Polyakov loop correlation 
functions calculated at finite temperature plays a simi¬ 
lar important role in the discussion of properties of the 
strong force as the static quark potential does at zero 
temperature. 

The properties of this observable (at T = 0: poten¬ 
tial, at T yf 0: free energy) at short and intermediate 
distances {rT < 1) is important for the understanding of 
in-medium modifications of heavy quark bound states. 
A quantitative analysis of heavy quark free energies be¬ 
comes of considerable importance for the discussion of 
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possible signals for the quark gluon plasma formation 
in heavy ion collision experiments |^|5|. For instance, 
recent studies of heavy quarkonium systems within po¬ 
tential models use the quark anti-quark free energy to 
define an appropriate finite temperature potential which 
is considered in the non-relativistic Schrodinger equation 
QHiQ Such calculations, however, do not quite 
match the results of direct lattice calculations of the 
quarkonium dissociation temperatures which have been 
obtained so far only for the pure gauge theory [13 im 
It was pointed out that the free energy (F) of a 
static quark anti-quark pair can be separated into two 
contributions, the internal energy (U) and the entropy 
(S'). The separation of the entropy contribution from 
the free energy, i.e. the variable U = F + TS, could 
define an approp riate effective potential at finite tem¬ 
perature^ 14ff(D 2^) = U, to be used as input in 

model calculations and might explain in parts the quanti¬ 
tative differences found when comparing solutions of the 
Schrodinger equation with direct calculations of spectral 
functions |l(illl| . First calculations which use the inter¬ 
nal energy obtained in our calculations [I1II1II3II3 
support this expectation. Most of these studies con¬ 
sider so far quenched QCD. Using potentials from the 
quenched theory, however, will describe the interaction 
of a heavy quark anti-quark pair in a thermal medium 
made up of gluons only. It is then important to under¬ 
stand how these results might change for the case of a 


^ While a definition of the quark anti-quark potential can be given 
properly at zero temperature using large Wilson loops, at finite 
temperature a definition of the thermal modification of an ap¬ 
propriate po tent ial energy between the quark anti-quark pair is 
complicated M- 
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thermal heat bath which also contains dynamical quarks. 

On the other hand, it is the large distance property of 
the heavy quark interaction which is important for our 
understanding of the bulk properties of the QCD plasma 
phase, e. q. the screening property of the quark gluon 
plasma |^, the equation of stat e 1^ and the 
order parameter (Polyakov loop) Elllfllll. In all 
of these studies deviations from perturbative calculations 
and the ideal gas behavior are expected and were indeed 
found at temperatures which are only moderately larger 
than the deconfinement temperature. This calls for quan¬ 
titative non-perturbative calculations. Also in this case 
most of todays discussions of the bulk thermodynamic 
properties of the QGP and its apparent deviations from 
the ideal gas behavior rely on results obtained in lat¬ 
tice studies of the pure gauge theory, although several 
qualitative differences are to be expected when taking 
into account the influence of dynamical fermions; for in¬ 
stance, the phase transition in full QCD will appear as 
an crossover rather than a ’true’ phase transition with re¬ 
lated singularities in thermodynamic observables. More¬ 
over, in contrast to a steadily increasing confinement 
interaction in the quenched QCD theory, in full QCD 
the strong interaction below deconfinement will show a 
qualitatively different behavior at large quark anti-quark 
separations. Due to the possibility of pair creation the 
stringlike interaction between the two test quarks can 
break leading to a constant potential and/or free energy 
already at temperatures below deconfinement |2(i| . 

Thus it is quite important to extend our recently de¬ 
veloped concepts for the analysis of the quark anti-quark 
free energies and internal energies in pure gauge theory 
[iM 1^ 1^ to the more complex case of QCD with 
dynamical quarks, and to quantify the qualitative differ¬ 
ences which will show up between pure gauge theories 
and QCD. 
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TABLE I: Sample sizes at each /3 value and the temperature 
in units of the (pseudo-) critical temperature Tc. 

For our study of the strong interaction in terms of 
the quark anti-quark free energies in full QCD lattice 
configurations were generated for 2-flavor QCD {Nf=2) 
on 16^ X 4 lattices with bare quark mass ma=0.1, i.e. 
m/T=0A, corresponding to a ratio of pion to rho masses 
(rriTr/mp) at the (pseudo-) critical temperature of about 




FIG. 1: (a) The color singlet quark anti-quark free energies, 
Fi{r, T), at several temperatures close to the phase transition 
as function of distance in physical units. Shown are results 
from lattice studies of 2-flavor QCD. The solid line represents 
in each figure the T — 0 heavy quark potential, V(r). The 
dashed error band corresponds to the string breaking energy 
at zero temperature, V (cbreaking) — 1000 — 1200 MeV, based 
on the estimate of the string breaking distance, rbreaking — 
1.2—1.4 fm [^ . (b) The color averaged free energy, Fnqjr, T), 
normalized such that Fav{r,T) = Fqq{r,T) — Tin9 ll^ ap¬ 
proaches the heavy quark potential, V (r) (line), at the small¬ 
est distance available on the lattice. The symbols are chosen 
as in (a). 


0.7 (a denotes the lattice spacing) [^. We have used 
Symanzik improved gauge and p4-improved staggered 
fermion actions. This combination of lattice actions is 
known to reduce the lattice cut-off effects in Polyakov 
loop correlation functions at small quark anti-quark sep¬ 
arations seen as an improved restoration of the broken 
rotational symmetry. For any further details of the sim¬ 
ulations with these actions see . In Table Q] we 

summarize our simulation parameters, i.e. the lattice 
coupling /3, the temperature T/T^ in units of the pseudo 
critical temperature and the number of configurations 
used at each /3-value. The pseudo critical coupling for 
this action is /3c = 3.649(2) [^. To set the physical 
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scale we use the string tension, cra^, measured in units 
of the lattice spacing, obtained from the large distance 
behavior of the heavy quark potential calculated from 
smeared Wilson loops at zero temperature |^. This is 
also used to define the temperature scale and a\fo is 
used for setting the scale for the free energies and the 
physical distances. For the conversion to physical units, 
^ = 420MeV is used. For instance, we get = 202(4) 
MeV calculated from Td\/a = 0.48(1) [^. In parts of 
our analysis of the quark anti-quark free energies we are 
also interested in the flavor and finite quark mass depen¬ 
dence. For this reason we also compare our 2-flavor QCD 
results to the todays available recent findings in quenched 
^=0) [III 23 and 3-flavor QCD [rn-nlmp ~ 0.4 
|34l| . Here we use Tc = 270 MeV for quenched and 
Tc = 193 MeV for the 3-flavor case. 

Our results for the color singlet quark anti-quark free 
energies, Fi, and color averaged free energies, Fq^, are 
summarized in Fig. function of distance at several 
temperatures close to the transition. At distances much 
smaller than the inverse temperature (rT 1) the dom¬ 
inant scale is set by distance and the QCD running cou¬ 
pling will be controlled by the distance. In this limit 
the thermal modification of the strong interaction will 
become negligible and the finite temperature free energy 
will be given by the zero temperature heavy quark poten¬ 
tial (solid line). With increasing quark anti-quark sep¬ 
aration, however, thermal effects will dominate the be¬ 
havior of the finite temperature free energies [rT 3> 1). 
Qualitative and quantitative differences between quark 
anti-quark free energy and internal energy will appear 
and clarify the important role of the entropy contribu¬ 
tion still present in free energies. The quark anti-quark 
internal energy will provide a different look on the inter¬ 
quark interaction and thermal modifications of the finite 
temperature quark anti-quark potential. Further details 
of these modifications on the quark anti-quark free and 
internal energies will be discussed. 

This paper is organized as follows: We start in sec¬ 
tion m with a discussion of the zero temperature heavy 
quark potential and the coupling. Both will be calcu¬ 
lated from 2-flavor lattice QCD simulations. We ana¬ 
lyze in section llTTI the thermal modifications on the quark 
anti-quark free energies and discuss quarkonium binding. 
Section llVl contains our summary and conclusions. A de¬ 
tailed discussion of the quark anti-quark internal energy 
and entropy will be given separately [s^. 


II. THE ZERO TEMPERATURE HEAVY 
QUARK POTENTIAL AND COUPLING 




FIG. 2: (a) The heavy quark potential at T = 0 from 1^ 

obtained from 2-flavor QCD lattice simulations with quark 
masses ma = 0.1 for different values of the lattice coupling (3. 
Fig. 2(b) shows an enlargement of the short distance distance 
regime. The data are matched to the bosonic string potential 
(dashed line) at large distances. Included is also the fit to the 
Cornell form (solid line) given in Eq. @. Note here that the 
heavy quark potential from quenched lattice QCD and the 
string model potential coincide already at r^/a > 0.8 [Mil 
(r > 0.4 fm). 

To eliminate the divergent self-energy contributions we 
matched these data for all /3-values (different /3-values 
correspond to different values of the lattice spacing a) at 
large distances to the bosonic string potential, 

V(r) = 


( 1 ) 


A. Heavy quark potential at T = 0 

For the determination of the heavy quark potential at 
zero temperature, V (r), we have used the measurements 
of large smeared Wilson loops given in [^ for the same 
simulation parameters {Nf=2 and ma = 0.1) and action. 


where we already have separated the Casimir factor so 
that Ostr = 7r/16. In this normalization any diver¬ 
gent contributions to the lattice potential are eliminated 
uniquely. In Fig. [21 we show our results together with the 
heavy quark potential from the string picture (dashed 
line). One can see that the data are well described by 
Eq. CJ at large distances, i.e. r^/a > 0.8, corresponding 
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to r > 0.4 fm. At these distances we see no major differ¬ 
ence between the 2-flavor QCD potential obtained from 
Wilson loops and the quenched QCD potential which can 
be well parameterized within the string model already for 
r > 0.4 fm 1^1^. In fact, we also do not see any signal 
for string breaking in the zero temperature QCD heavy 
quark potential. This is expected due to the fact that the 
Wilson loop operator used here for the calculation of the 
T = 0 potential has onl y sm all overlap with states where 
string breaking occurs Moreover, the distances 

for which we analyze the data for the QCD potential are 
below r < 1.2 fm at which string breaking is expected to 
set in at zero temperature and similar quark masses |29| . 


B. The coupling at T = 0 

Deviations from the string model and from the pure 
gauge potential, however, are clearly expected to be¬ 
come apparent in the 2-flavor QCD potential at small 
distances and may already be seen from the short dis¬ 
tance part in Fig. [3 These deviations are expected to 
arise from an asymptotic weakening of the QCD cou¬ 
pling, i.e. a = a{r), and to some extent also due to the 
effect of including dynamical quarks, i.e. from leading 
order perturbation theory one expects 

87r/3olog(l/(rAQCD)) ’ 


with 


/3o 


33 - 2Nf 
487r2 ’ 


(3) 


where Nf is the number of flavors and Aqcd denotes the 
corresponding QCD-A-scale. The data In Fig. Hb) show 
a slightly steeper slope at distances below ~ 0.5 
compared to the pure gauge potential given in Ref. 
indicating that the QCD coupling gets stronger in the 
entire distance range analyzed here when Including dy¬ 
namical quarks. This Is in qualitative agreement with 
©■ To include the effect of a stronger Coulombic part 
in the QCD potential we test the Cornell parameteriza¬ 
tion, 


V{r) 4: a 

s/a 3 r/a 


(4) 


with a free parameter a. From a best-fit analysis of 
Eq. 0) to the data ranging from 0.2 < r/a < 2.6 we 
find 


a = 0.212(3) . (5) 

This already may indicate that the logarithmic weaken¬ 
ing of the coupling with decreasing distance will not too 
strongly Influence the properties of the QCD potential at 
these distances, i.e. at r > 0.1 fm. However, the value 
of a Is moderately larger than Ogtr — 0.196 Introduced 


above. To compare the relative size of a in full QCD 
to a In the quenched theory we again have performed a 
best-fit analysis of the quenched zero temperature poten¬ 
tial given in using the Ansatz given in Eq. 0 and a 
similar distance range. Here we find Oquenched = 0.195(1) 
which is again smaller than the value for the QCD cou¬ 
pling but quite comparable to Ostr- In earlier studies of 
the heavy quark potentials In pure gauge theories and full 
QCD even larger values for the couplings were reported 
IS Ql- To avoid here any confusions 

concerning the value of a we should stress that a should 
not be mixed with some value for the QCD coupling con¬ 
stant aqcD, it simply Is a fit parameter indicating the 
’average strength’ of the Coulomb part in the Cornell po¬ 
tential. The QCD coupling could be identified properly 
only In the entire perturbative distance regime and will 
be a running coupling, i.e. aqcD = Q!qcd(?')- 

When approaching the short distance perturbative 
regime the Cornell form will overestimate the value of 
the coupling due to the perturbative logarithmic weak¬ 
ening of the latter, aqcD = aqcD(?')- To analyze the 
short distance properties of the QCD potential and the 
coupling in more detail, i.e. for r < 0.4 fm, and to 
firmly establish here the onset of its perturbative weak¬ 
ening with decreasing distance, it is customary to do so 
using non-perturbative definitions of running couplings. 
Follo wing the discussions on the running of the QCD cou¬ 
pling |4^ |4^ 1^ [ 43 , it appears most convenient to 
study the QCD force, i.e. dV{r)/dr, rather than the 
QCD potential. In this case one defines the QCD cou¬ 
pling in the so-called qq-scheme. 


agq{r) 


3^2 dV{r) 
4 dr 


( 6 ) 


In this scheme any undetermined constant contribution 
to the heavy quark potential cancels out. Moreover, the 
large distance, non-perturbative confinement contribu¬ 
tion to aqq{r) is positive and allows for a smooth match¬ 
ing of the perturbative short distance coupling to the 
non-perturbative large distance confinement signal. In 
any case, however. In the non-perturbative regime the 
value of the coupling will depend on the observable used 
for Its definition. 

We have calculated the derivatives of the potential with 
respect to the distance, dV{r)/dr, by using finite differ¬ 
ence approximations for neighboring distances on the lat¬ 
tice for each /3-value separately. Our results for aqq (r) as 
a function of distance in physical units for 2-flavor QCD 
are summarized In Fig. |21 The symbols for the /3-values 
are chosen as in Fig. laa). We again show in that fig¬ 
ure the corresponding line for the Cornell fit (solid line). 
At large distances, r > 0.4 fm, the data clearly mimic 
the non-perturbative confinement part of the QCD force, 
aqq{r) ~ Sr^afA. We also compare our data to the recent 
high statistics calculation in pure gauge theory (thick 
solid line) [^. These data are available for r >0.1 fm 
and within the statistics of the QCD data no significant 
differences could be Identified between the QCD and pure 
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FIG. 3: The short distance part of the running coupling 
aqq ir) in 2-flavor QCD at zero temperature defined in Eq. @ 
as function of the distance r (in physical units). The symbols 
for the different /3-values are chosen as indicated in Fig. Ha). 
The lines are discussed in the text. 

gauge data for r > 0.4 fm. At smaller distances (r < 0.4 
fm), however, the data show some enhancement com¬ 
pared to the coupling in quenched QCD. The data below 
0.1 fm, moreover, fall below the large distance Cornell 
fit. This may indicate the logarithmic weakening of the 
coupling. At smaller distances than 0.1 fm we therefore 
expect the QCD potential to be influenced by the weak¬ 
ening of the coupling and aqq{r) will approach values 
clearly smaller than a deduced from the Cornell Ansatz. 
Unfortunately we can, at present, not go to smaller dis¬ 
tances to clearly demonstrate this behavior with our data 
in 2-flavor QCD. Moreover, at small distances cut-off ef¬ 
fects may also influence our analysis of the coupling and 
more detailed studies are required here. Despite these 
uncertainties, however, in earlier studies of the coupling 
in pure gauge theory US m Ei it is shown that the 
perturbative logarithmic weakening becomes already im¬ 
portant at distances smaller than 0.2 fm and contact with 
perturbation theory could be established. 

As most of our lattice data for the finite tempera¬ 
ture quark anti-quark free energies do not reach distances 
smaller than 0.1 fm we use in the following the Cornell 
form deduced in 0 as reference to the zero temperature 
heavy quark potential. 


III. QUARK ANTI-QUARK FREE ENERGY 

We will analyze here the temperature dependence of 
the change in free energy due to the presence of a heavy 
(static) quark anti-quark pair in a 2-flavor QCD heat 
bath. The static quark sources are described by the 
Polyakov loop, 

(7) 


with 

W{x) = Y[Uo{x,t) , ( 8 ) 

r —1 

where we already have used the lattice formulation with 
Uo{x, t) e SU{3) being defined on the lattice link in time 
direction. The change in free energy due to the presence 
of the static color sources in color singlet (Ai) and color 
octet (Fs) states can be calculated in terms of Polyakov 
loop correlation functions HElliail, 

g-F,(T)/T+C ^ iTr(W(x)lUt(y)) , (9) 

O 

g-FsW/T-HC ^ l{TTW{x)TTW\y)) - 

8 

lTT{Wix)W\y)) , (10) 

where r = \x — y\. As it stands, the correlation functions 
for the color singlet and octet free energies are gauge 
dependent quantities and thus gauge fixing is needed to 
define them properly. Here, we follow |49l| and fix to 
Coulomb gauge. In parts we also consider the so-called 
color averaged free energy defined through the manifestly 
gauge independent correlation function of two Polyakov 
loops, 

f.-Fqq(r)/T+c ^ i(TrtU(f)TrtUt(0)) 

9 

= {L{x)L\y)) . (11) 

The constant C appearing in ®, (Cni) and cu also 
includes divergent self-energy contributions which re¬ 
quire renormalization. Following |12| the free energies 
have been normalized such that the color singlet free en¬ 
ergy approaches the heavy quark potential (solid line) at 
the smallest distance available on the lattice, Fi{r/a = 
1,T) = V{r). In Sec. IIII5I we will explain the connec¬ 
tion of this procedure to the the renormalized Polyakov 
loop and show the resulting renormalization constants in 
Table CH 

Some results for the color singlet, octet and averaged 
quark anti-quark free energies are shown in Fig. H for 
one temperature below and one temperature above de¬ 
confinement, respectively. The free energies calculated 
in different color channels coincide at large distances and 
clearly show the effects from string breaking below and 
color screening above deconfinement. The octet free en¬ 
ergies above Tc are repulsive for all distances while below 
Tq the distances analyzed here are not small enough to 
show the (perturbatively) expected repulsive short dis¬ 
tance part. Similar results are obtained at all temper¬ 
atures analyzed here. In the remainder of this section 
we study in detail the thermal modifications of these free 
energies from short to large distances. We begin our anal¬ 
ysis of the free energies at small distances in Sec. nun 
with a discussion of the running coupling which leads 
to the renormalization of the free energies in Sec. mill 


L(f) = iTrVF(f) 

o 
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FIG. 4: Heavy quark free energies for 2 flavors of dynamical 
quarks at a quark mass of m/T = 0.4 calculated on 16^ x 4 
lattices. Shown are the free energies in different color chan¬ 
nels, the singlet (Fi), octet (Fs) and color averaged {Fqq) 
free energies normalized here as discussed in m to the zero 
temperature potential obtained in Sec. im (solid line). 



FIG. 5: The running coupling in the gg-scheme defined in 
Eq. 1141 calculated from derivatives of the color singlet free 
energies with respect to r at several temperatures as function 
of distance below and above deconfinement. We also show the 
corresponding coupling at zero temperature (solid line) from 
Eq. 0 and compare the results again to the results in pure 
gauge theory (thick solid and dashed lines) 


The separation of small and large distances which char¬ 
acterizes sudden qualitative changes in the free energy 
will be discussed in Sec. imn Large distance modifica¬ 
tions of the quark anti-quark free energy will be studied 
in Sec. inmi at temperatures above and in Sec. CHI at 
temperatures below deconfinement. 

Our analysis of thermal modifications of the strong in¬ 
teraction will mainly be performed for the color singlet 
free energy. In this case a rather simple Coulombic r- 
dependence is suggested by perturbation theory at T = 0 
and short distances as well as for large distances at high 
temperatures. In particular, a proper r-dependence of 
Fqq is difficult to establish . This maybe is attributed 
to contributions from higher excited states or to the 
repulsive contributions from states with static charges 
fixed in an octet configuration. 


A. The running coupling at T 7^ 0 


We extend here our studies of the coupling at zero 
temperature to finite temperatures below and above de¬ 
confinement following the conceptual approach given in 
[^. In this case the appropriate observable is the color 
singlet quark anti-quark free energy and its derivative. 
We use the perturbative short and large distance rela¬ 
tion from one gluon exchange , i. e. in the limit 

?'Aqcd ^ 1 zero temperature perturbation theory sug¬ 
gests 


4 cy(r] 

F,{r,T) ^ V{r) ^ -3-^, (12) 


while high temperature perturbation theory, i. e. rT ^ 1 
and T well above Tc, yields 

Fi(r,T) ~ . (13) 

O V 

In both relations we have neglected any constant contri¬ 
butions to the free energies which, in particular, at high 
temperatures will dominate the large distance behavior of 
the free energies. Moreover, we already anticipated here 
the running of the couplings with the expected dominant 
scales r and T in both limits. At finite temperature we 
define the running coupling in analogy to T = 0 as (see 

in [23), 




3^^ dF,{r,T) 
4 dr 


(14) 


With this definition any undetermined constant con¬ 
tributions to the free energies are eliminated and the 
coupling defined here at finite temperature will re¬ 
cover the coupling at zero temperature defined in 
in the limit of small distances. Therefore aqq{r,T) will 
show the (zero temperature) weakening in the short dis¬ 
tance perturbative regime. In the large distance limit, 
however, the coupling will be dominated by Eq. ca 
and will be suppressed by color screening, aqq{r,T) ~ 
a(T) exp(— mD(T)r), rT ^ 1. It thus will exhibit a 
maximum at some intermediate distance. Although in 
the large distance regime aqq{r,T) will be suppressed 
by color screening and thus non-perturbative effects will 
strongly control the value of aqq{r,T), in this limit the 
temperature dependence of the coupling, a{T), can be 
extracted by directly comparing the singlet free energy 
with the high temperature perturbative relation above 
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deconfinement. Results from such an analysis will be 
given in Sec. Em 

We calculated the derivative, dFi/dr, of the color sin¬ 
glet free energies with respect to distance by using cu¬ 
bic spline approximations of the r-dependence of the free 
energies for each temperature. We then performed the 
derivatives on basis of these splines. Our results for 
aqq{r,T) calculated in this way are shown in Fig. E and 
are compared to the coupling at zero temperature dis¬ 
cussed already in Sec. nrm Here the thin solid line cor¬ 
responds to the coupling in the Cornell Ansatz deduced 
in Eq. 0. We again show in this hgure the results 
from SU (3)-lattice (thick line) and perturbative (dashed 
line) calculations at zero temperature from 1^ . The 
strong r-dependence of the running coup ling near Tc ob¬ 
served already in pure gauge theory is also visible in 
2-flavor QCD. Although our data for 2-flavor QCD do not 
allow for a detailed quantitative analysis of the running 
coupling at smaller distances, the qualitative behavior is 
in quite good agreement with the recent quenched results. 
At large distances the running coupling shows a strong 
temperature dependence which sets in at shorter sepa¬ 
rations with increasing temperature. At temperatures 
close but above aqq{r,T) coincides with aqq{r) al¬ 
ready at separations r ~ 0.4 fm and clearly mimics here 
the confinement part of aqq{r). This is also apparent in 
quenched QCD [^. Remnants of the confinement part 
of the QCD force may survive the deconfinement transi¬ 
tion and could play an important role for the discussion of 
non-perturbative aspects of quark anti- quar k interactions 
at temperatures moderately above Tc [ 13 . A clear 

separation of the different effects usually described by the 
concept of color screening (T > Tc) and effects usually 
described by the concept of string-breaking (T < Tc) is 
difficult to establish at temperatures in the close vicinity 
of the confinement deconfinement crossover. 

We also analyzed the size of the maximum that the 
running coupling aqq{r,T) at fixed temperature exhibits 
at a certain distance, rmax, *-e. we identify a temperature 
dependent coupling, aqq{T), defined as 



FIG. 6: The size of the maximum, &qq(T), defined in Eq. II51 . 
as function of temperature in 2-flavor QCD (fille d symbols) 
and pure gauge theory (open symbols) from 1^ . The lines 
are explained in the text. 


and quenched QCD become visible in aqq{T). Nonethe¬ 
less, the value of aqq{T) drops from about 0.5 at temper¬ 
atures only moderately larger than the transition tem¬ 
perature, T > 1.2Tc, to a value of about 0.3 at 2Tc. This 
change in aqq{T) with temperature calculated in 2-flavor 
QCD does not appear to be too dramatic and might in¬ 
deed be described by the 2-loop perturbative coupling. 




(16) 


with 


Po 

Pi 


1 


1 


11 - 


(IOtt^)^ 




The values for r^ax will be discussed in Sec ITTT71 (see 
Fig.E). Values for diqq{T) are also available in pure gauge 
theory at temperatures above deconfinement^. Our 
results for aqq{T) in 2-flavor QCD and pure gauge the¬ 
ory are shown in Fig.E^s function of temperature, T /Tc. 
At temperatures above deconfinement we cannot identify 
signihcant differences between the data from pure gauge 
and 2-flavor QCD^. Only at temperatures quite close but 
above the phase transition small differences between full 


^ In pure gauge theory rmax and aqq(T) would be infinite below 
Tc. 

® Note here, however, the change in temperature scale from Tc = 
202 MeV in full and Tc = 270 MeV in quenched QCD. 


assuming vanishing quark masses. In view of the am¬ 
biguity in setting the scale in perturbation theory, /tT, 
we performed a best-fit analysis to fix the scale for 
the entire temperature range, 1.2 < T/Tc < 2. We 
find here p, = 1.14(2)^ with Tc/A^g^ = 0.77(21) u sing 
Tc ~ 202(4) MeV lH and ~ 261(17) MeV [53||, 
which is still in agreement with the lower limit of the 
range of scales one commonly uses to fix perturbative 
couplings, = TT, ...,Att. This is shown by the solid line 
(fit) in Fig.Eincluding the error band estimated through 
/i = TT to ^ = 47r and the error on Tc/Aj^ (dotted 
lines). We will turn back to a discussion of the tempera¬ 
ture dependence of the coupling above deconfinement in 
Sec. mTPl 

At temperatures in the vicinity and below the phase 
transition temperature, T < 1.2Tc, the behavior of 
UqqpT) is, however, quite different from the perturba¬ 
tive logarithmic change with temperature. The values 


















for aqq{T) rapidly grow here with decreasing tempera¬ 
ture and approach non-perturbatively large values. This 
again shows that aqq{r,T) mimics the confinement part 
of the zero temperature force still at relatively large dis¬ 
tances and that this behavior persists up to temperatures 
close but above deconfinement. This again demonstrates 
the persistence of confinement forces a.t T > Tc and in¬ 
termediate distances and demonstrates the difficulty to 
separate clearly the different effects usually described by 
color screening and string breaking in the vicinity of the 
phase transition. We note here, however, that similar to 
the coupling in quenched QCD the coupling which 
describes the short distance Coulombic part in the free 
energies is almost temperature independent in this tem¬ 
perature regime, i.e. even at relatively large distances 
the free energies shown in Fig. ^ show no or only little 
temperature dependence below deconfinement. 

B. Renormalization of the quark anti-quark free 
energies and Polyakov loop 


breaking (T < Tc). These plateau values, Foo{T) = 
Fi{r oo,T), are decreasing with increasing temper¬ 
ature in the temperature range analyzed here. In general 
it is expected that F'oo(T) will continue to increase at 
smaller temperature and will smoothly match V(r = oo) 
at zero temperature while it will become negative at 
high temperature and asymptotically is expected to be¬ 
come proportional to g^T 1^ . The plateau yalue of 
the quark anti-quark free energy at large distances can 
be used to define non-perturbatively the renormalized 
Polyakov loop ^3) *-6- 

(18) 

As the unrenormalized free energies approach |(L)p at 
large distances, this may be reinterpreted in terms of a 
renormalization constant that has been determined by 
demanding to hold at short distances 

L-- ^ \{{Zig,m))^^ L)\ . (19) 


On the lattice the expectation value of the Polyakov 
loop and its correlation functions suffer from linear di¬ 
vergences. This leads to vanishing expectation values in 
the continuum limit, a —> 0, at all temperatures. To 
become a meaningful physical obseryable a proper renor¬ 
malization is required [l3. miHi. We follow here the 
conceptual approach suggested in [H M and extend 
our earlier studies in pure gauge theory to the present 
case of 2-fiayor QCD. First experiences with this renor¬ 
malization method in full QCD were already reported in 

[2II2I. 

In the limit of short distances, r <C 1/T, thermal modi¬ 
fications of the quark anti-quark free energy become neg¬ 
ligible and the running coupling is controlled by distance 
only. Thus we can fix the free energies at small distances 
to the heayy quark potential, Fi{r 1/T,T) ~ Vir), 
and the renormalization group equation (RGE) will lead 
to 


limT 

r— 


dF^{r,T) 

dT 


0 , 


(17) 


where we already haye assumed that the continuum limit, 
a ^ 0, has been taken. On basis of the analysis of the 
coupling in Sec. mm and experiences with the quark 
anti-quark free energy in pure gauge theory [l2.1^ we 
assume here that the color singlet free energies in 2-flayor 
QCD calculated on finite lattices with temporal extent 
Nr = 4 already haye approached appropriate small dis¬ 
tances, r <C 1/T, allowing for renormalization. 

The (renormalized) color singlet quark anti-quark free 
energies, Fi{r,T), and the heayy quark potential, V{r) 
(line), were already shown in Fig. ^a) as function of 
distance at seyeral temperatures close to the phase tran¬ 
sition. From that figure it can be seen that the quark 
anti-quark free energy fixed at small distances approaches 
finite, temperature dependent plateau yalues at large dis¬ 
tances signaling color screening (T > Tc) and string 


The yalues for Z{g, m) for our simulation parameters are 
summarized in Tabled The normalization constants for 
the free energies appearing in mni are then giyen by 

C = -2NrZ{g,m). (20) 

An analysis of the renormalized Polyakov loop expecta¬ 
tion value in high temperature perturbation theory 
suggests at (resummed) leading order^, the behavior 

L-"(r) l + ^!^^a(T) (21) 

in the fundamental representation. Thus high tem¬ 
perature perturbation theory suggests that the limiting 
value at infinite temperature, L’'®°(T —> 00 ) = 1 is ap¬ 
proached from above. An expansion of m then suggests 
Foo(T) ~ —^mjj(T)a{T) ~ —O(g^T). We thus expect 
Foc(T) —>■ —00 in the high temperature limit. 

To avoid here any fit to the complicated r- and T- 
dependence of the quark anti-quark free energy we esti¬ 
mate the value of Foc{T) from the quark anti-quark free 
energies at the largest separation available on a finite 
lattice, r = Nrj/2. As the free energies in this renormal¬ 
ization scheme coincide at large distances in the different 
color channels we determine F,^(T) from the color av¬ 
eraged free energies, i.e. Fco{T) = Fqq(r = N„/2,T). 
This is a manifestly gauge invariant quantity. In Fig. [7| 


In Ref. the Polyakov loop expectation value is calcu¬ 

lated in pure gauge theory and the Debye mass, mrtiT)/T = 
\/Nc/3g(T), enters here through the resummation of the gluon 
polarization tensor. When changing from pure gauge to full QCD 
quark loops will contribute to the polarization tensor. In this case 
resummation will lead to the Debye mass given in II24I . Thus the 
flavor dependence in Eq. 121 1 at this level is given only by the 
Debye mass. 
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p 

Z{g,m) 

TjTc 

^ren(y) 

3.52 

1.333(19) 

0.76 

0.033(2) 

3.55 

1.351(10) 

0.81 

0.049(2) 

3.60 

1.370(08) 

0.90 

0.093(2) 

3.63 

1.376(07) 

0.96 

0.160(3) 

3.65 

1.376(07) 

1.00 

0.241(5) 

3.66 

1.375(06) 

1.02 

0.290(5) 

3.68 

1.370(06) 

1.07 

0.398(7) 

3.72 

1.374(02) 

1.16 

0.514(3) 

3.75 

1.379(02) 

1.23 

0.575(2) 

3.80 

1.386(01) 

1.36 

0.656(2) 

3.85 

1.390(01) 

1.50 

0.722(2) 

3.90 

1.394(01) 

1.65 

0.779(1) 

3.95 

1.396(13) 

1.81 

0.828(3) 

4.00 

1.397(01) 

1.98 

0.874(1) 

4.43 

1.378(01) 

4.01 

1.108(2) 


TABLE II: Renormalization constants, Z{g,m), versns /3 and 
the renormalized Polyakov loop, versus T/Tc for 2-flavor 
QCD with quark mass m/T = 0.4. 


we show the results for L''®" in 2-flavor QCD (filled sym¬ 
bols) compared to the quenched results (open symbols) 
obtained in In quenched QCD is zero below 
Tc as the quark anti-quark free energy signals permanent 
confinement, i.e. FcoiT < Tc) = oo in the infinite volume 
limit, while it jumps to a finite value just above Tc. The 
singularity in the temperature dependence of lo¬ 

cated at Tc clearly signals the first order phase transition 
in SU(3) gauge theory. The renormalized Polyakov loop 
in 2-flavor QCD, however, is no longer zero below Tc- 
Due to string breaking the quark anti-quark free ener¬ 
gies approach constant values leading to non-zero values 
of T’'®". Although the renormalized Polyakov loop cal¬ 
culated in full QCD is no longer an order parameter for 
the confinement deconfinement phase transition, it still 
shows a quite different behavior in the two phases and 
a clear signal for a qualitative change in the vicinity of 
the transition. Above deconfinement T"'®"(T) yields finite 
values also in quenched QCD. 

In the temperature range 1 < T/Tc < 2 we find that 
in 2-flavor QCD T“'®" lies below the results in quenched 
QCD. This, however, may change at higher tempera¬ 
tures. The value for T’’®'' at 4Tc is larger than unity and 
we find indication for > TqTenched(47’c)- 

The properties of T''®'', however, clearly depend on the 
relative normalization of the quark anti-quark free ener¬ 
gies in quenched and full QCD. 


C. Short vs. large distances 

Having discussed the quark anti-quark free energies at 
quite small distances where no or only little temperature 
effects influence the behavior of the free energies and at 



FIG. 7: The renormalized Polyakov loop in 2-flavor QCD 
(filled circles) compared to the quenched results (open sym¬ 
bols) from 1121 as function of temperature in units of the (dif¬ 
ferent) transition temperatures. 

quite large distances where aside from T no other scale 
controls the free energy, we now turn to a discussion of 
medium effects at intermediate distances. The aim is 
to gain insight into distance scales that can be used to 
quantify at which distances temperature effects in the 
quark anti-quark free energies set in and may influence 
the in-medium properties of heavy quark bound states in 
the quark gluon plasma. 

It can be seen from Fig. Hla) that the color singlet free 
energy changes rapidly from the Coulomb-like short dis¬ 
tance behavior to an almost constant value at large dis¬ 
tances. This change reflects the in-medium properties of 
the heavy quark anti-quark pair, i.e. the string-breaking 
property and color screening. To characterize this rapid 
onset of in-medium modifications in the free energies we 
introduced in Ref. 0 a scale, rmedi defined as the dis¬ 
tance at which the value of the T = 0 potential reaches 
the value Too(T), i.e. 

V{r^cd) = Too(T) . (22) 

As T’oo(T) is a gauge invariant observable this relation 
provides a non-perturbative, gauge invariant definition of 
the scale rmed- While in pure gauge theory the color sin¬ 
glet free energies signal permanent confinement at tem¬ 
peratures below Tc leading to a proper definition of this 
scale only above deconfinement, in full QCD it can be 
deduced in the whole temperature range. On the other 
hand, the change in the coupling aqq(r,T) as function 
of distance at fixed temperature mimics the qualitative 
change in the interaction when going from small to large 
distances and the coupling exhibits a maximum at some 
intermediate distance. The location of this maximum, 
fmax, can also be used to identify a scale that charac¬ 
terizes separation between the short distance vacuumlike 
and the large distance medium modified interaction be¬ 
tween the static quarks Due to the rapid crossover 
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FIG. 8 : The location of the maximum in aqq{r,T) at 

fixed temperature, Vmax {Nf—O: open squares, Nf=2-. filled 
squares), and our results for {Nf=0'. open circles, Nf=2: 
filled circles, Nf=3'. crosses) defined in Eq. 1221 versus T/Tc- 
The band at the left frame indicates the distance range at 
which string breaking is expected to occur in 2-flavor QCD at 
T = 0 and quark mass m-Klvrip ~ 0.7 [^ . The various lines 
are explained in the text. 


from short to large distance behavior (see Fig. da)) it 
should be obvious that r^ed and r^ax define similar 
scales, however, by construction r^nax Tmed- 

To gain important information about the flavor and 
quark mass dependence of our analysis of the scales in 
QCD, we also took data for Fr^jT) from Ref. at 
smaller quark mass, m-Tr/mp ~ 0.4 |^, and calculated 
Tmed in 3-flavor QCD with respect to the parameteriza¬ 
tion of V{r) given in |^. It is interesting to note here 
that a study of the flavor and quark mass dependence 
of Cmed and rmax is independent of any undetermined 
and maybe flavor and/or quark mass dependent overall 
normalization of the corresponding V (r) at zero tempera¬ 
ture. Our results for Vmax (-^/=0,2) and r^ed {Nf=0,2,3) 
are summarized in Fig.|Hlas function of T/Tc. It can be 
seen that the value r^nax — 0.6 fm is approached almost 
in common in quenched and 2-flavor QCD at the phase 
transition and it commonly drops to about 0.25 fm at 
temperatures about 2Tc. No or only little differences be¬ 
tween r^nax calculated from pure gauge and 2-flavor QCD 
could be identified at temperatures above deconfinement. 
The temperature dependence of r^ed is similar to that of 
i^max and again we see no major differences between pure 
gauge {Nf=0) and QCD {Nf=2,3) results. In the vicin¬ 
ity of the transition temperature and above both scales 
almost coincide. In fact, above deconfinement the flavor 
and finite quark mass dependence of rmed appears quite 
neg ligible. At high temperature we expect rmed — ^/gT 
Il2l| while in terms of Vmax we found agreement with 
rmax = 0.48(1) fm Tr /T (solid lines) at temperatures 
ranging up to 12Tc [^. Note that both scales clearly lie 
well above the smallest distance attainable by us on the 
lattice, rT = l/N^ = 1/4. This distance is shown by the 


lower dashed line in Fig.|51 

At temperatures below deconfinement rmax and rmed 
rapidly increase and fall apart when going to smaller tem¬ 
peratures. In fact, at temperatures below deconhnement 
we clearly see difference between rmed calculated in 2 - 
and 3-flavor QCD. To some extend this is expected due 
to the smaller quark mass used in the 3-flavor QCD study 
as the string breaking energy gets reduced. It is, how¬ 
ever, difficult to clearly separate here a finite quark mass 
effect from flavor dependence. In both cases rmed ap¬ 
proaches, already at T ~ 0.8Tc, quite similar values to 
those reported for the distance where string breaking at 
T = 0 is expected at similar quark masses. In 2-flavor 
QCD at T = 0 and quark mass TOtt/top ~ 0.7 the string 
is expected to break at about 1.2 — 1.4 fm while at 
smaller quark mass, uiT^Imp ~ 0.4 it might break earlier 

M- 

In contrast to the complicated r- and T-dependence 
of the free energy at intermediate distances high tem¬ 
perature perturbation theory suggests a color screened 
Coulomb behavior for the singlet free energy at large dis¬ 
tances. To analyze this in more detail we show in Fig. El 
the subtracted free energies, r(Fi(oo,T) — Fi{r,T)). It 
can be seen that this quantity indeed decays exponen¬ 
tially at large distances, rT > 1. This allows us to study 
the temperature dependence of the parameters a(T) and 
moiT) given in Eq. 11311 . At intermediate and small 
distances, however, deviations from this behavior are ex¬ 
pected and can clearly be seen and are to some extent 
due to the onset of the r-dependence of the coupling at 
small distances. These deviations from the simple ex¬ 
ponential decay become important already below some 
characteristic scale, r^, which we can roughly identify 
here as r^T ~ 0.8 — 1. This scale which defines a lower 
limit for the applicability of high temperature perturba¬ 
tion theory is shown by the upper dashed line in Fig. |H| 
(r^T = 1). It lies well above the scales rmed and rmax 
which characterize the onset of medium modifications on 
the quark anti-quark free energy. 


D. Screening properties above deconfinement and 
the coupling 

1. Screening properties and quarkonium binding 

We follow here the approach commonly used pcl 
miii and define the non-perturbative screening mass, 
niDir), and the temperature dependent coupling, a(T), 
from the exponential fall-off of the color singlet free en¬ 
ergies at large distances, rT > 0.8 - 1. A consistent 
definition of screening masses, however, is accompanied 
by a proper definition of the temperature dependent cou¬ 
pling and only at sufficiently high tem per atures contact 
with perturbation theory is expected . A simi¬ 

lar discussion of the color avera ged quark anti-quark free 
energy is given in Refs, [liil^l^. 

We used the Ansatz 03 to perform a best-fit analysis 
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FIG. 9: The color singlet free energy versus rT obtained from 
lattice calculations in 2-flavor QCD at several temperatures 
above deconfinement. 


of the large distance part of the color singlet free energies, 
i.e. we used fit functions with the Ansatz 

Ai(r,T) - Fi(r = oo,r) = (23) 

3r 

where the two parameters a{T) and m{T) are used to 
estimate the coupling a(T) and the Debye mass mniT), 
respectively. The fit-range was chosen with respect to 
our discussion in Sec. imn i.e. rT > 0.8 - 1, where we 
varied the lower fit limit within this range and averaged 
over the resulting values. The temperature dependent 
coupling a(T) defined here will be discussed later. Our 
results for the screening mass, mu{T)/T, are summa¬ 
rized in Fig. ^las function of T/Tc and are compared to 
the results obtained in pure gauge theory . The data 
obtained from our 2-flavor QCD calculations are some¬ 
what larger than in quenched QCD. Although we are not 
expecting perturbation theory to hold at these small tem¬ 
peratures, this enhancement is in qualitative agreement 
with leading order perturbation theory, i.e. 


mpiT) 

T 


N f 


1/2 


9{T). 


(24) 


However, using the 2-loop formula to estimate the 
temperature dependence of the coupling leads to signif¬ 
icantly smaller values for mp/T even when setting the 
scale by /t = TT which commonly is used as an upper 
bound for the perturbative coupling. We therefore fol¬ 
low MM and introduce a multiplicative constant. A, 
i.e. we allow for a non-perturbative correction defined as 


mp{T) 

T 


= A 



ff2 — loop(T) , 


(25) 


and fix this constants by best agreement with the non- 
perturbative data for mp{T)/T at temperatures T > 1.2. 
Here the scale in the perturbative coupling is fixed by 


FIG. 10: The non-perturbatively defined screening masses 
from the large distance behavior of Fi{r,T) calculated using 
16® X 4 lattices as function of temperature in 2-fiavor QCD 
(filled circles) and in pure gauge theory (open squares Im I 
using 32® X 4 lattices. The solid lines show the fit given in 
Eq. I23II with the corresponding error band (dotted lines). 


/i = 27r. This analysis leads to A = 1.417(19) and is 
shown as solid line with error band (dotted lines). Simi¬ 
lar results were already reported in |l9t 1^ for screening 
masses in pure gauge theory. Using the same fit range, 
i.e. T = 1.2Tc — 4.1Tc, for the quenched results, we ob¬ 
tain A = 1.515(17). To avoid here any confusion concern¬ 
ing A we note that its value will crucially depend on the 
temperature range used to determine it. When approach¬ 
ing the perturbative high temperature limit, A —> 1 is 
expected. 

It is interesting to note here that the difference in 
mp/T apparent in Fig. between 2-flavor QCD and 
pure gauge theory disappears when converting mp{T) 
to physical units. This is obvious from Fig. ^2 which 
shows the Debye screening radius, rp = limp. In gen¬ 
eral rp is used to characterize the distance at which 
medium modifications of the quark anti-quark interac¬ 
tion become dominant. It often is used to describe the 
screening effects in phenomenological inter-quark poten¬ 
tials at high temperatures. From perturbation theory one 
expects that the screening radius will drop like IjgT. A 
definition of a screening radius, however, will again de¬ 
pend on the ambiguities present in the non-perturbative 
definition of a screening mass, mp{T). A different quan¬ 
tity that characterizes the onset of medium effects, rmed, 
has already been introduced in Sec. irrm this quantity 
is also expected to drop like l/{gT) at high tempera¬ 
tures and could be considered to give an upper limit for 
the screening radius [^. In Fig. we compare both 
length scales as function of temperature, T/Tc, and com¬ 
pare them to the findings in quenched QCD mill 
It can be seen that in the temperature range analyzed 
here rp(T) < rmed{T) and no or only little differ¬ 
ences between the results from quenched {Nf=Q) and full 
{Nf=2,3) QCD could be identified. Again we stress that 
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FIG. 11: The screening radius estimated from the inverse 
Debye mass, vd = llmo {Nf=0-. open squares, Nf=2 filled 
squares), and the scale r^ed open circles, Nf=2: 

filled circles, A^/=3: crosses) defined in 1221 as function of 
T/Tc- The dotted line indicates the smallest separation avail¬ 
able on lattices with temporal extent Nt = 4. The horizontal 
lines gives the mean squared charge radii of some charmonium 
states, J/tp, Xc and ip! (see also [l3l) and the band at the left 
frame shows the distance at which string breaking is expected 
in 2-fiavor QCD at T = 0 and quark mass ra-^jrap ~ 0.7 1^ . 


in the perturbative high temperature limit differences are 
expected to arise as expressed by Eq. CT . 

It is important to realize that at distances well be¬ 
low r^ed medium effects become suppressed and the 
color singlet free energy almost coincides with the zero 
temperature heavy quark potential (see Fig. da)). In 
particular, the screening radius estimated from the in¬ 
verse Debye mass corresponds to distances which are only 
moderately larger than the smallest distance available in 
our calculations (compare with the lower dotted line in 
Fig. IHJ). In view of the almost temperature independent 
behavior of the color singlet free energies at small dis¬ 
tances (Fig. da)) it could be misleading to quantify the 
dominant screening length of the medium in terms of 
ty) = llmu- On the other hand the color averaged free 
energies show already strong temperature dependence at 
distances similar to ro (see Fig.db)). 

Following we also included in Fig. ^2 the mean 
charge radii of the most prominent charmonium states, 
Xc and ipp as horizontal lines. These lines char¬ 
acterize the averaged separation r which enters the ef¬ 
fective potential in potential model calculations. It thus 
is reasonable to expect that the temperature at which 
these radii equal r^ed could give a rough estimate for 
the onset of thermal effects in the charmonium states. 
It appears quite reasonable from this view that J/pj in¬ 
deed may survive the phase transition Eim, while Xc 
and ipt are supposed to show significant thermal modifi¬ 
cations at temperatures close to the transition. Recent 
potential model calculations support this analysis |l(i| . 
The wave functions for these states, however, will also 



FIG. 12: The temperature dependence of the coupling a{T) 
as function of temperature (filled circles) from fits of the large 
distance behavior of free energies using Eq. We also 

show again the values of aqq (T) defined as the maximum value 
of the running coupling aqq (r, T) (filled squares) discussed in 
Sec. midl and the perturbative 2-loop coupling, Eq. GUi with 
scales /r = TT,..., dvr (dashed lines). The open symbols indicate 
results from corresponding quenched QCD calculations 1^ . 

reach out to larger distances and this estimate can 
only be taken as a first indication for the relevant tem¬ 
peratures. Further details on this issue including also 
bottomonium states have been given in Ref. 0 We 
will turn again to a discussion of thermal modifications 
of quarkonium states in Ref. using finite temperature 
quark anti-quark energies. 


2. Temperature dependence of Qs 

We finally discuss here the temperature dependence of 
the QCD coupling, a{T), extracted from the fits used 
to determine also mo, i-e. from Eq. ll^ . From fits of 
the free energies above deconfinement we find the values 
shown in Fig. E|as function of T/Tc given by the filled 
circles. We again show in this figure also the tempera¬ 
ture dependent coupling aqq{T) introduced in Sec. lIIIfO 
It can clearly be seen that the values for both couplings 
are quite different, aqq{T) > a(T), at temperatures close 
but above deconfinement while this difference rapidly de¬ 
creases with increasing temperature. This again demon¬ 
strates the ambiguity in defining the coupling in the 
non-perturbative temperature range due to the different 
non-perturbative contributions to the observable used for 
its definition [^. In fact, at temperatures close to the 
phase transition temperature we find quite large values 
for a{T), i.e. a{T) ~ 2 — 3 in the vicinity of Tc, while it 
drops rapidly to values smaller than unity, i.e. a(T) < 1 
already at temperatures T/T^ > 1.5. A similar behav¬ 
ior was also found in for the coupling in pure gauge 
theory (open symbols). In fact, no or only a marginal 
enhancement of the values calculated in full QCD com- 
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FIG. 13: The plateau value of the quark anti-quark free en¬ 
ergy, FooiT), calculated in 2-flavor QCD as function of T/Tc 
at temperatures in the vicinity and below the phase transi¬ 
tion. The dashed lines show the expected value at zero tem¬ 
perature, F(rbreaking) , with Tbreaking 1.2 ~ 1.4 fm USiug 
quark mass m-^lrrip ~ 0.7 1^ . The open symbols (T > Tc) 
correspond to Foa{T) in quenched QCD (from Ref. [Ty]) and 
the crosses to 3-fiavor QCD studies (from 1^ 1. The relative 
normalization of the corresponding V (r) used for renormal¬ 
ization of Fi{r —> 00,T) in quenched and full QCD is such 
that the Cornell parameterization of V (r) does not contain 
any constant at large distances. 


pared to the values in quenched QCD could be identified 
here at temperatures T < 1.5Tc. We stress again that the 
large values for a{T) found here should not be confused 
with the coupling that characterizes the short distance 
Coulomb part of Fi{r,T). The latter is almost temper¬ 
ature independent at small distances and can to some 
extent be described by the zero temperature coupling. 


E. String breaking below deconfinement 


can be calculated from the lightest heavy-light meson, 
^'breaking — 1.2 — 1.4 fm and is shown on the left side 
in Fig.^within the dotted band. It can be seen that r^ed 
in 2-flavor QCD does indeed approach such values at tem¬ 
peratures T < 0.8Tc. This suggests that the dependence 
on temperature in 2-flavor QCD is small below the small¬ 
est temperature analyzed here, 0.76Tc. This can also be 
seen from the behavior of Foo{T) shown in Fig. (see 
also Fig. ma)) compared to the value commonly expected 
at T = 0. We use 17(rbreaking) — 1000— 1200 MeV as ref¬ 
erence to the zero temperature string breaking energy 
with quark mass m-^lmp ~ 0.7. This estimate is shown 
on the left side in Fig. El as the dotted band. A similar 
behavior is expected for the free energies in 3-flavor QCD 
and smaller quark mass, m^/rup ~ 0.4. As seen also in 
Fig.Elthe values for Fao{T) are smaller than in 2-flavor 
QCD and larger quark mass. This may indicate that 
string breaking sets in at smaller distances for smaller 
quark masses. However, in [30) no mass dependence (in 
the color averaged free energies) was observed below the 
quark mass analyzed by us (m/T=0.4). At present it is, 
however, difficult to judge whether the differences seen 
for 2- and 3-flavor QCD for T/Tc < 1 are due to quark 
mass or flavor dependence of the string breaking. Al¬ 
though Foc{T) still is close to 17(rbreaking) at T ^ 0.8Tc, 
it rapidly drops to about half of this value in the vicin¬ 
ity of the phase transition, Fco{Tc) ~ 575 MeV. This 
value is almost the same in 2- and 3-flavor QCD; we And 
F^^='^{Tc) ~ 575(15) MeV and F^f=^{Tc) ~ 548(20) 
MeV. It is interesting to note that also the values of 
Foo{T) in quenched QCD [Nf=0) approach a similar 
value at temperatures just above Tc- We And Fao{T ^) ~ 
481(4) MeV where T+ = 1.02rc denotes the closest tem¬ 
perature above Tc analyzed in quenched QCD. Of course, 
the value for Fco {T ^) will increase when going to tem¬ 
peratures even closer to Tc- The flavor and quark mass 
dependence of FooiT) including also higher temperatures 
will be discussed in more detail in Ref. |^. 


We finally discuss the large distance properties of the 
free energies below Tc. In contrast to the quark anti¬ 
quark free energy in quenched QCD where the string be¬ 
tween the quark anti-quark pair cannot break and the 
free energies are linearly rising at large separations, in 
full QCD the string between two static color charges can 
break due to the possibility of spontaneously generat¬ 
ing gg-pairs from the vacuum. Therefore the quark anti¬ 
quark free energy reaches a constant value also below Tc. 
In Fig. ^this behavior is clearly seen. 

The distances at which the quark anti-quark free ener¬ 
gies approach an almost constant value move to smaller 
separations at higher temperatures. This can also be 
seen from the temperature dependence of Tmed in Fig. H 
at temperature below Tc. By construction r^ed describes 
a distance which can be used to estimate a lower limit 
for the distance where the string breaking will set in. 
An estimate of the string breaking radius at T = 0 


IV. SUMMARY AND CONCLUSIONS 

Our analysis of the zero temperature heavy quark po¬ 
tential, V(r), calculated in 2-flavor lattice QCD using 
large Wilson loops shows no signal for string break¬ 
ing at distances below 1.3 fm. This is quite consis¬ 
tent with earlier findings |2^ The r-dependence 

of V{r) becomes comparable to the potential from the 
bosonic string picture already at distances larger than 
0.4 fm. Similar findings have also been reported in lat¬ 
tice studies of the potential in quenched QCD . 

At those distances, 0.4 fm < r < 1.5 fm, we find no 
or only little differences between lattice data for the po¬ 
tential in quenched (V/=0) given in Ref. and full 
iNf=2) QCD. At smaller distances, however, deviations 
from the large distance Coulomb term predicted by the 
string picture, Ostr — 0.196, are found here when per¬ 
forming best fit analysis with a free Cornell Ansatz. We 
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find a ~ 0.212(3) which could describe the data down to 
r > 0.1 fm. By analyzing the coupling in the gg-scheme 
defined through the force, dV{r)/dr, small enhancement 
compared to the coupling in quenched QCD is found for 
r < 0.4 fm. At distances substantially smaller than 0.1 
fm the logarithmic weakening of the coupling enters and 
will dominate the r-dependence of V{r). The observed 
running of the coupling may already signal the onset of 
the short distance perturbative regime. This is also evi¬ 
dent from quenched QCD lattice studies of V (r) . 

The running coupling at finite temperature defined in 
the gg-scheme using the derivative of the color singlet 
quark anti-quark free energy, dFi{r,T)/dr, shows only 
little qualitative and quantitative differences when chang¬ 
ing from pure gauge to full QCD at temperatures 
well above deconfinement. Again, at small distances the 
running coupling is controlled by distance and becomes 
comparable to cxqq{r) at zero temperature. The prop¬ 
erties of aqq(r,T) at temperatures in the vicinity of the 
phase transition are to large extent controlled by the con¬ 
finement signal at zero temperature. A clear separation 
of the different effects usually described by the concepts 
of color screening {T > Tc) and the concept of string 
breaking (T < Tc) is difficult in the crossover region. 
Remnants of the confinement part of the QCD forces 
may in parts dominate the non-perturbative properties of 
the QCD plasma at temperatures only moderately larger 
than Tc- This supports similar findings in recent studies 
of the quark anti-quark free energies in quenched QCD 

isg. 

The properties of the quark anti-quark free energy 
and the coupling at small distances thus again allow for 
non-perturbative renormalization of the free energy and 
Polyakov loop 0. The crossover from confinement to 
deconfinement is clearly signaled by the Polyakov loop 
through a rapid increase at temperatures close to Tc- 
String breaking dominates the quark anti-quark free en¬ 
ergies at temperatures well below deconfinement in all 
color channels leading to finite values of the Polyakov 
loop. The string breaking energy, Foo(T), and the dis¬ 
tance where string breaking sets in, are decreasing with 
increasing temperatures. The plateau value Fao{T) ap¬ 
proaches about 95% of the value one usuall y es timates 
at zero temperature, P(r'breaking) — 1-1 GeV al¬ 

ready for T ~ O.STc- We thus expect that the change in 
quark anti-quark free energies is only small when going 
to smaller temperatures and the quark anti-quark free 
energy, T'i(r, T), will show only small differences from 
the heavy quark potential at T = 0, V{r). Significant 
thermal modifications on heavy quark bound states can 
thus be expected only for temperatures above O.STc. Our 
analysis of r^ed suggests indeed a qualitative similar be¬ 
havior for the free energies in 3-flavor QCD. This can also 
be seen from the behavior of Tmed shown in Fig. 1111 

At temperatures well above the (pseudo-) critical tem¬ 
perature, i.e. 1.2 < T/Tc < 4, no or only little 
qualitative differences in the thermal properties of the 
quark anti-quark free energies calculated in quenched 


{Nf=0) and full {Nf=2,3) QCD could be established 
here when converting the observables to physical units. 
Color screening clearly dominates the quark anti-quark 
free energy at large distances and screening masses, which 
are non-perturbatively determined from the exponential 
fall-off of the color singlet free energies, could be ex¬ 
tracted (for Nf=2). In accordance with earlier findings 
in quenched QCD [T^l^ we find substantially larger val¬ 
ues for the screening masses than given by leading order 
perturbation theory. The values of the screening masses, 
m£)(T), again show only marginal differences as function 
oiTjTc compared to the values found in quenched QCD 
(see also Fig. EH). The large screening mass defines a 
rather small screening radius, td = which refers 

to a length scale where the singlet free energy shows al¬ 
most no deviations from the heavy quark potential at 
zero temperature. It thus might be misleading to quan¬ 
tify the length scale of the QCD plasma where temper¬ 
ature effects dominate thermal modifications on heavy 
quark bound states with the observable td = IIvtid in 
the non-perturbative temperature regime close but above 
Tc- On the other hand the color averaged free energies 
show indeed strong temperature dependence at distances 
which could be characterized by \jmB- In view of color 
changing processes as a mechanism for direct quarkonium 
dissociation the discussion of the color averaged free 
energy could become important. 

We have also compared ru and Tmed in Fig. the 
expected mean squared charge radii of some charmonium 
states. It is reasonable that the temperatures at which 
these radii equal Tmed give a first indication of the tem¬ 
perature at which thermal modihcations become impor¬ 
tant in the charmonium states. It appears thus quite rea¬ 
sonable that J/'i/' will survive the transition while Xc and 
il)i are expected to show strong thermal effects at tem¬ 
peratures in the vicinity of the transition and this may 
support recent findings [Ulllil. Of course the wave 
functions of these states will also reach out to larger dis¬ 
tances and thus our analysis can only be taken as a first 
indication of the relevant temperatures. We will turn 
back to this issue in Ref. |^. The analysis of bound 
states using, for instance, the Schrodinger equation will 
do better in this respect. It can, however, clearly be 
seen from Fig.^]that although rmed{Tc) — 0.7 fm is ap¬ 
proached almost in common for A^y=0,2,3, it falls apart 
for Nf=2,3 at smaller temperatures. It thus could be dif¬ 
ficult to determine suppression patterns from free ener¬ 
gies for quarkonium states which are substantially larger 
than 0.7 fm independently from Nj and/or hnite quark 
mass. 

The analysis presented here has been performed for a 
single quark mass value that corresponds to a pion mass 
of about 770 MeV (m^/mp ~ 0.7). In Ref. js^, how¬ 
ever, no major quark mass effects were visible in color 
averaged free energies below this quark mass value. The 
comparisons of rmed and Foo{T) calculated in 2-flavor 
{niTr/mp ~ 0.7) with results calculated in 3-fiavor QCD 
(m^/TOp ~ 0.4 [^) supports this property. While at 
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temperatures above deconfinement no or only little dif¬ 
ferences in this observable can be identified, at tempera¬ 
tures below Tc differences can be seen. To what extend 
these are due to the smaller quark masses used in the 
3-flavor case or whether these differences reflect a fla¬ 
vor dependence of the string breaking distance requires 
further investigation. The present analysis was carried 
out on one lattice size (16^ x 4) and therefore perform¬ 
ing an extrapolation to the continuum limit could not be 
done with the current data. However the analysis of the 
quenched free energies ,where no major differences 

between the Nt = 4 and Nr = 8 results were visible, and 
the use of improved actions suggests that cut-off effects 
might be small. Despite these uncertainties and the fact 
that parts of our comparisons to results from quenched 
QCD are on a qualitative level, we find quite important 
information for the study of heavy quark bound states in 
the QCD plasma phase. At temperatures well above 
i.e. 1.2 < T/Tc < 4, no or only little differences appear 
between results calculated in quenched and QCD. This 
might suggest that using thermal parameters extracted 
from free or internal energy in quenched QCD as input for 
model calculations of heavy quark bound states [iMlS i® 
a reasonable approximation. Furthermore this also sup¬ 
ports the investigation of heavy quarkonia in quenched 
lattice QCD calculations using the analysis of spectral 
functions On the other hand, however, most 


of our 2- and 3-flavor QCD results differ from quenched 
calculations at temperatures in the vicinity and below 
the phase transition. Due to these qualitative differences, 
results from quenched QCD could make a discussion of 
possible signals for the quark gluon plasma production in 
heavy ion collision experiments complicated when tem¬ 
peratures and/or densities close to the transition become 
important. 
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